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INTRODJCTIjN 


Tha  note  submitted  her-.-./ith  deals  with  the  control  of  the 
directivity  pattern  of  sound  radiatec"  from  a  plate  vibrating  in  an 
infinite  ba.ffle.   It  supplements  recent  investigations  of  the  pro- 
blem such  as  thosJ  by  the  Harvard  Underwater  Sound  Laboratory  and 
by  H.  I'yquisty  in  so  far  as  complete  mathematical  solution  of  the 
following  pro'-lem  is  given!  How  can  a.  prescribed  radiation  pattern 
be  obtained  by  a  suitable  distribution  of  amplitudes  (and  phases) 
ovc^r  a  vibrating  plate?  The  general  answer,  obtained  for  cir- 
cular plates  in  a.n  explicit  form,  c:.n,  of  course,  be  realized  in 
practice  only  to  a  more  or  less  approximate  degree.   Nevertheless, 
the  idealized  solution  given  herein  will  serve  as  an  indice.tion 
of  how  the  designing  engineer  could  profitably  distribute  sia- 
plitudes  (and  phases)  over  the  various  elements  of  a  mosaic  of 
vibrating  elements  v/hich  togcth'-r  form  the  v/hole  of  a  projector, 

R.  Coiirant 

Official  Invest ieat or 


•5^  "Report  on  Directivity  Pattern  of  Sound  Sources"  O.S.R.D.  No.  7O6 
"Directivity  of  Soueed  in  Water  -  Elementary  Arrays"  C4  -  MDRC  -  O72 


RADIATION  PATTBRIJ3  FROM  A  VIBRATING  PLATE 

by  E,  Bromberg  and  H,  Cohn 

The  problem  of  the  control  of  the  directivity  pattern  of  sound 
radiated  from  a  vibrator  set  in  an  infinite  rigid  wall  or  "baffle"  has 
recently*  been  treated  by  the  Harvard  Underwater  Sound  Laboratory  and 
by  Dr.  Harry  Nyquist  of  Division  C  Section  4,   Their  vibrators  con- 
sisted of  line  sources  and  point  sources,** 

Tfe  wish  to  consider  the  same  problem  for  a  two-dimensional 
array  of  sound  sources  so  denseljr  packed  that  it  may  be  considered 
as  a  continuous  surface  distribution  of  a  vibrating  plate.  Vfe  assume 
further  that  the  overall  dimensions  of  the  plate  and  the  frequency 
of  vibration  have  been  specified. 

We  take  our  plate  in  the  x-y  plane  vibrating  in  the  '^  z- 
direction,  so  that  the  eqxjation  of  the  surface  of  the  plate  at  any 
time  t  is 
(1)  z  =  Re  ':  f(x,y)  e  '^"^^^ 

where  f(x,y)  is  a  given  amplitude  function,  extending  over  the 
projection  of  the  plate  on  the  x-y  plane.  If  all  points  on  the 
plate  vibrate  in  phase  (or  out  of  phase  by  180°)  f(x,y)   is  real 
and  positive  (or  negative);  otherwise  f(x,y)   is  complex.  We  shall 
study  the  dependence  of  the  radiation  pattern  on  the  amplitude  func- 
tion f(x,y). 


*  "Report  on  Directivity  Pattern  of  Sound  Sources"  O.S.R.D.  No.  706 

"Directivity  of  Sound  in  mter  -  Elementary  Arrays"  04-  -  N.D.R.C.  -  072 
+*  Further  investigation  of  the  regulation  of  line  sources  is  being 
made  by  Mr.  R.  P.  Shaiv  of  the  N.  Y.  U.  Applied  Mathematics  Group, 
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The   complex  amplitude   function  of   the  acoustic  pressure,   at 
the   point  of  observation,   P(x      ,   y   ,    Zq),    is   given  hy: 


(2) 


"■CiUkr 


P     -    ^£^1   I  &_f(x,y)dx  dy 


Zrr- 


fr 


plate 


where 


-  V  / 


\/   (x-x^)'^  +■  iy-y^)''  ^-(Zq)' 


In  what  follovi^s  we    confine   our  attention  to   a   circular  plate  where 


f(x,y)  ^-f(R), 


2  ^  „Z 


R  =■  V  X  -^'  y 


0<R<  A,. 


where  A  is  the  radius  of  the  plate. 
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Wc  now  consider  formula  (a)  at  large  distances  from  the  plate.   Then 


r  ^;:.r  ■\-   y  sin.'^ 


where  r   -  \l  ^n     +"  ^o^  '+'^0^  *  '^o  ^^  *^°  angle  beti/veen  the  radial 
vector  to  (xq,  y^,  z^)  and  the  z  axis,  and  the  x,y  axes  are  chosen, 
for  convenience,  so  that  x  —  0,   Hence   (2)   may  be  simplified: 
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e       J     J     e  f  (  V  x2  -H  y2)    dx  dy 


o 

where  wg  define     g     "by 


;^)     -   J^J   J      e''^^f(V^2T7F)    dx   dy 


x^  -t-  y  V  A^ 


^^  .-27r 


^  f(H)   R  dR    \      e  '^  ^  Cltt' 


2lT 
/.A 


(3)    '^    g^)     =      \   R  f(R)   Jq  (/A.R)  dR 
R-6 


ikr^ 


The   complex  amplitude  of  the  wave  is   thengiven  "by  z^lf^^jiTi/Vj—-: const, 

.1  " 

and  its   intensity  "by     const.  X  |  g(  ^^  sin^J  f       The   "amplitude  radiation 

pattern"  is   defined  by 

/O^)^    IgC^sin^)! 

If  the  plate  is  rigid,  that  is,   f(R)  =  constant.^Tj^,  v;e  have 
the  well  known  result 

(4)  g(  ,CC)  ^  VI  a2  ll^Vi 

'^  A/U. 

Let  us  now  suppose  thr.t  the  radiation  pattern  is  preassigned, 
i.e.,  the  non-negative  function  y{^)      is  given  for  0</f7<2,«   T^-is 
defines   g(A^   for  0<y-(<K(g  may  be  taken  as  zero  elsewhere).   Then 
if  we  find  a  "mode  of  vibration"  f(R)   such  that 
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(3a)  n 

g{jL)    '   \  rJ^(^^r)  -f  (r)  dn 
o 

this  implies  that  an  infinite  vibrating  plate  with  amplitude  function 
f(R)  would  produce  the  desired  radiation  pattern  V^(,i9^) ,  'Ye  may  in- 
Tert  (.3a)  by  the  Fourier-Bcssel  Theorem,  v/hich  states  that 

(5)  g(yx)  =  J  R  J^(P./-)  dRJf  s(|)  j;(Rf)  df 

R^O  ^--"^ 

Hence,   from  (3a),   since  g(  ^  )   ="0     for     ^  >  ^    we  set 

k 

(6)  -f(R)-   J      S(|)  J      (f  R)fc('f 

f-'O 

for  all  R  from  0  tooo. 

The  f amotion  given  in  (6)  scorns  to  represent  a  vibrating  plate 
of  infinite  radius,  but  in  practice  \ve  shall  use  a  pl?vte  of  finite 
radius  A.   Hence  instead  of  achieving  the  desired  radiation  pattern  of 
(3a),  we  shall  find  actually. 


(3b) 


gA  (/^)  -     /  '^J^il^^)ii^)   dR 


The  difference  betyreen  (3a),  the  desired  radiation  pattern,  and 
(3b),  the  actual  radiation  pattern,  becomes  negligible  for  large  plates 
by  the  convergence  of  the  integral  in  (3a).   The  closeness  of  the 
approximation  to  a  desired  radiation  pattern,  achieved  by  setting  the 
amplitude  of  vibration  equal  to  f (R)  only  over  the  disk  of  radius  A, 
will  then  be  limited  by  the  permissible  magnitudes  of  A  and  k, 
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EXAJvIPLES 


Suppose  the  following  radiation  pattern  is  desired; 

^yO^)   -  1      for      0  <^<^ 
('  /{^)=    0      for     4;^-'-^^? 


Cs  (/-) "  1 

for              0  </-^    <Mo 

(   ^{/i.)   -0 

for             /(^  </^ 

f(R)  = 

that  is. 


Then,  from  (6) 


and  the  v/ay  in  which  this  function  applies  to  the  preceding  discussion 
is  obvious . 

--^-      I    2h 
b)      Since  the  Fourier  Bessel  Transform  of  e  '^    is  -—  fc   ,  they 

h 

will  he  proportional  to  g( /Lc)  and  f(R)  respectiTely.   It  is  clear  that 
if  f(R)  is  sensibly  different  fron  zero  only  over  a  very  small  range 
of  R,  the  resulting  beam  will  be  very  vn.de;  on  the  other  hand  the 
beam  can  be  sharpened  to  any  desired  extent  by  increasing  the  range 
of  R  over  which  f(R)  is  appreciably  different  from  zero  —  which  would 
involve  increasing  the  radius  of  the  disk  accordingly.   The  "reci- 
procity" of  this  pattern  has  been  noticed  in  the  report  previously 
referred  to  as  O.3.R.D.  No.  706 
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APPLI CATIONS 


In  most  practical  cases  the  construction  of  the  vibrating  plate 
is  not  such  as  to  permit  any  arbitrary  continuous  variation  of  the 
complex  amplitude  of  vibration  over  the  surface  of  the  plate. 

For  the  disk,  the  simplest  case  is  that  in  v/hich  the  complex 
amplitude  is  held  constant  over  the  entire  surface.   The  "amplitude 
radiation  pattern"  of  the  resulting  beam  was  expressed  in  (4)  as 

f.^                     rL    .     0^^                 ?  J.CAksinO^) 
(4)       g(ksin.t7)  -  TT-^  A'^    _! 

A  i<  s  in  .6^ 

For  purposes  of  computation  vre   have  taken  the  radius  of  the  disk 
A-  STTcm,  k  ^^  1  per  cm,  and  we  adjust  v^  so  that  the  maximum  amplitude 
of  g  is  unity.  The  resulting  function  has  been  plotted  on  the  ac- 
companying chart  as  curve  a.   If  we  consider  as  negligible  any  parts 
of  the  beam  which  have  an  amplitude  of  less  than  one-tenth  of  the 
maximum,  vre   see  that  this  vibrating  disk  yields  a  bean  with  a  strong 
central  conical  portion  together  with  "side-lobes"  vAich  extend  ap- 
proximately to  an  angle  of  39'-'  with  the  axis  of  the  disk.  In  v/hat 
follows  ;ve  shall  take  as  our  criterion  of  a  desirable  change,  one  which 
will  decrease  the  angle  of  the  principal  conical  portion,  and  which 
will  elimiriate  the  "side-lobes".   It  is  assumed  that  the  radius  of 
the  disk  and  K  are  to  remain  unchanged. 

It  has  been  found  possible  in  practice  to  break  the  disk  up 
into  a  central  circular  portion  and  concentric  annuli,  with  the  com- 
plex amplitude  of  vibration  constant  over  each  section,  but  differing 
from  one  section  to  the  next.   This  involves,  in  effect,  an  attempt 
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to  develop  an  approximation  to  a  desirable  "mode  of  vibration"  function 

f(R). 

7fe  consider,  now,  the  case  of  a  disk  with  a  central  circle  of 
radius  B  and  a  single  annulus  of  outer  radius  A.   If  we  now  sot 
""v,    for     R<  B 
f  (R)  = 


I  1 
J 


V    for  B  <  R  <  A 
the  integration  of  (3)  yields 


(7)    g(  k  sin^)  =  V,  a2  J.  (-  ^  ^^^^^^  ,  (v,-v-)b2  JS^^-'-r.^) 

A  K  sin^9>  Bk  sina- 

t  V  a2  /''J^^i^'^^*^  ^-  (II  -  1)  B^j;(Bksin^) 

^    I   '^Aksln";^     "^2     I^    BK  sin -6^ 


The  parameter  v„  is  to  be  adjusted  so  that  the  maximum  amplitude  of  the 
beam  remains  unity,  as  in  the  previous  case.   This  leaves  us  with  the 
parameters  v^  and  B  by  means  of  which  we  can  control  the  directional 
properties  of  the  beam. 

As  an  example,  v/e  have  evaluated  g  for  the  case  B  -  0.546  A, 
and  V,  =  2.67  v  .   The  result  is  plotted  on  the  chart  as  curve  b_.   It 
can  be  seen  that  the  "side-lobes"  have  been  almost  completely  elimina- 
ted, but  that  the  central  portion  has  become  somewhat  wider  than  in 
curve  a.   Further  consideration  indicates  that  the  further  narrowing 
of  the  central  portion  of  curve  b  v/ill  lead  to  an  increase  in  the 
amplitude  of  the  "side-lobes"  and  that  the  central  beam  can  never  be 
narrower  than  that  of  curve  a. 

If  the  disk  were  to  be  divided  into  three  or  four  concentric 
sections,  it  seems  reasonable  to  expect  that  the  beam  could  be 
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improved  even  more,  because  of  the  increase  in  the  namber  of  paramet- 
ers. Hov/ever,  it  is  not  readily  apparent  just  hov:  to  vary  these  para- 
meters to  achieve  the  best  results,  and  therefore  a  good  deal  but  not 
an  excessive  amount  of  computational  \vork  might  be  involved  in  this 
procedure.   Thus  the  expression  for  a  disk  with  tvo  annuli,  '.vith  C 
the  innermost  radius,  B  the  intermediate  radius  and  A  the  outermost 
radius,  and  with 


(-1 

) 
< 


for     R  <  C 


f(R)-  <  V  for  C<P  <  B 
3  for  B  <  R  <  A 


leads   to   the   following  x'alue   for   g: 

2   "^1^^^  sine)   vg  B^   H^k  sine) 

(8)   g(Ksin.S^)  zr  VgA     Ak  sinC      "^ %  '  ^  ^^  "Tk^Ii^ 

/"^l  ■^2  1       C^        J^^CCk  sin3>) 

Again  v   is  adjusted  to  make  the  maximum  amplitude  of  the  beam  unity, 
leaving  v-,  ,  Vo,  C,  and  B  as  parameters  to  determine  the  directions.l 
properties  of  the  beam.  We  have  made  a  rough  attempt  to  improve  the 
beam  in  this  case,  keeping  the  v.  real,  however.   The  result  is  plotted 
as  curve  o  on  the  chart,  with  the  values  of  the  parameters  indicated 
there.  As  can  be  seen,  it  is  an  improvement  over  both  of  the  other 
cutnres.   Presumably  a  better  result  could  be  achieved  with  more  care, 
but  it  serves  to  indicate  the  advantages  of  this  method. 
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The  same  method  is  equally  applicable  to  vibrating  plates  with 
other  shapes,  such  as  rectangular  plates.   In  many  oases  it  is  possible 
to  divide  the  plate  up  into  a   "mosaic"  of  small  sections,  each  vi- 
brating with  its  own  complex  amplitude,  so  adjusted  as  to  yield  a  beam 
with  the  particular  qualities  which  may  be  desired.   This  sort  of  problem 
will  yield  to  the  kind  of  analysis  indicated  above.   The  smaller  the 
sections,  the  closer  it  is  possible  to  approximate  to  a  desirable 
"mode  of  vibration"  f(R),  and  hence  the  closer  the  approximation  to 
the  desired  radiational  beam. 

It  is  obvious  that  this  method  would  yield  considerably  better 
approximations  if  the  specified  size  and  frequency  of  vibration  of  the 
plate  are  increased. 
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Data   for  Accompanying   Chart; 


/ 


A- 


0. 


A^  5fr   Crr\. 


A ^ 

Curve  a;     For   ona-soction  disc  or   disc  with  constant  amplitude  of  vibra- 
tion. 


Curve  b:      For   two- section  disc: 
B    r  0.546     A 


\-'    '-''     ^2 


Curve  c:   For  three- section  disc: 
C  =  0.6   A 
B  =  0.8 
T^  =    0.82  V3 
Vg  -  0.22  V3 


t-^^ 
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